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Using an effective-field-theory (nonlinear a model) description of interacting electrons in a dis- 
ordered metal ring enclosing magnetic flux, we calculate the moments of the persistent current 
distribution, in terms of interacting Goldstone modes (diffusons and cooperons). At the lowest or 
Gaussian order we reproduce well-known results for the average current and its variance that were 
originally obtained using diagrammatic perturbation theory. At this level of approximation the cur- 
rent distribution can be shown to be strictly Gaussian. The nonlinear a model provides a systematic 
way of calculating higher-order contributions to the current moments. An explicit calculation for 
the average current of the first term beyond Gaussian order shows that it is small compared to the 
Gaussian result; an order-of-magnitude estimation indicates that the same is true for all higher-order 
contributions to the average current and its variance. We therefore conclude that the experimentally 
observed magnitude of persistent currents cannot be explained in terms of interacting diffusons and 
cooperons. 



I. INTRODUCTION 

The unexplained large magnitude of observed persis- 
tent currents in disordered gold and copper rings - 
with a size on the order of micrometers and at tempera- 
tures in the milli-Kelvin range — has been an outstand- 
ing problem in condensed matter physics over the past 
half decade. The phenomenon was already predicted by 
Hund.EI Due to the high degree of difficulty of the exper- 
iments there is only a very limited number of observa- 
tions that have been published n u In the case of clean 
GaAs rings, where the electrons move ballisticaJly, the 
experimento is in good agreement with theoryB When 
impurities are present however, the electrons move dif- 
fusively; in this case there is a discrepancy of .one or 
two orders of magnitude between the experimentstru and 
the theoretical predictions, especially in the case of the 
single-ring experiment u 

On the theoretical side there has recently been a lot of 
activity.LTEll However, in spite of all efforts, there is still 
no satisfactory and widely accepted explanation for the 
observed current magnitude. It seems that a resolution of 
the problem will necessarily involve taking into account 
the (screened) Coulomb interaction between electrons. 
This belief is based on the results of exact diagonalization 
studies for tight-binding models, showing that electron- 
electron interaction can significantly enhance the magni- 
tude of the current. The size of the systems that can 
be studied in this manner is however far to small for 
any firm conclusion to be drawn from it. Larger system 
sizes can be achieved by numerically solving the problem 
in Hartree-Fock approximation, a method that has re- 
cently become popular .□ This approach however neglects 
potentially important electronic correlations, and should 
therefore be treated with some suspicion. 

Analytic approaches are mostly based-on Green's func- 
tion techniques for disordered systems £3 a nice summary 



of the work performed along these lines can be found kx 
Rcf. ^3|. Two widely accepted results for the averageO 
and typicalEJ current are based on these methods, and 
involve the summation of an infinite series of ladder di- 
agrams. These resummations give rise to so-called dif- 
fusons, related to the classical diffusion mode associated 
with the conserved electron density, and cooperons, the 
time reversed equivalent of diffusons. Unfortunately, the 
results obtained along this route do not explain the ex- 
periments. At higher orders in perturbation theory it 
becomes very difficult to determine which classes of di- 
agrams are important_aad there has been some contro- 
versy on this subject JHjEj It is clear that a resolution of 
the question whether or not a theory in terms of inter- 
acting diffusons and cooperons is capable of explaining 
the experimental observations, requires an unambiguous 
method for generalizing the known results to higher or- 
ders in perturbation theory. 

A consistent theoretical description of the interplay 
between disorder and electron-electron interaction is a 
notoriously hard problem in solid state physics. In the 
context of localization, a mapping of the quantum many- 
body problem to a classical field theory — a matrix 
nonlinear a model that provides a systematic descrip- 
tion in terms of the slow Goldstone modes associated 
with a continuous symmetry of the order parameter — 
has proved very useful. The method was pioneered by 
Wegner,a and extended to, include electron-electron in- 
teractions by Finkel'stein.cJ A general and self-contained 
derivation csua be found in a recent review by Belitz and 
KirkpatrickEa 

The nonlinear a model provides the starting point for 
the theory presented in this paper. Here we will only 
give an outline of the derivation, but we will provide full 
details when appropriate, in particular since there are 
necessarily some technical differences with Ref. |2^. By 
analogy with the localization problem we expect that the 
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nonlinear a model will allow us to perform a perturbative 
analysis in a much more efficient manner than is possible 
using standard perturbation theory. The results however 
will be equivalent, and in principle one should be able to 
identify the corresponding classes of Feynman diagrams. 
This has been done explicitly in some cases for the local- 
ization problem. 

The outline of the paper is as follows. In Se c. [ll| we 
formulate the problem. Subsequently , in Sec. Ill 



we 



introduce the nonlinear a model and present a system- 
atic method for calculating moments of the current dis- 
tribution. In Sec. IV we evaluate all moments at the 



lowest or Gaussian order, and show that at this level of 
approximation the current distribution is strictly Gaus- 
sian. In Sec. [v| the first contribution to the average cur- 
rent beyond Gaussian order is calculated in full detail; 
an estimate is obtained for the magnitude of all higher- 
order contributions to the first and second moment of 
the current distribution; furthermore it is argued that 
a low-frequency divergence occurs at all orders beyond 
Gaussian. The paper ends with a discussion in Sec. VI 



where ts = L/vp is the time it takes an electron to 
travel around the ring at the Fermi velocity vp. When 
disorder is present the electron motion is diffusive, with 
an associated time scale to = L 2 /D, where D = vpljd 
is the diffusion coefficient and d — 3 the dimensionality. 
Thus the diffusive current scale is 
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Three experiments have been performed on disordered 
metal rings: (i) a multiring experiments in which the av- 
erage current per ring in an array of 10 7 rings was found 
to be on the order of Ip>] (h) a single- ring experimentcl 
in which much larger currents on the order of Ib were 
observed; (iii) a 30-ring experimentE that interpolates 
between the first two experiments. The second experi- 
ment in particular poses a challenge, since theory so far 
predicts a value for the typical current that is one or two 
orders of magnitude smaller than what has been observed 
experimentally. 



II. FORMULATION OF THE PROBLEM 



III. NONLINEAR a MODEL APPROACH 



Let us describe in detail the problem we consider. We 
will use units in which Ti = c = 1. The circumference of 
the ring is denoted by L and its transverse width by L±. 
When the ring is placed in a stationary magnetic field, 
an equilibrium current ,B 



(1) 



will flow in response along the circumference of the ring. 
Here Q is the free energy and 
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is the magnetic flux enclosed by the ring. The density of 
impurities is specified in terms of the mean free path £, 
which is large compared to the Fermi wavelength: kp£ ~ 
100. The width L± ~ I of the ring is much smaller than 
its circumference L; typically we have L/l ~ 100. Thus it 
is reasonable to assume that we can describe the behavior 
of the electrons in terms of slow diffusive modes. 

It is well known that in strictly one-dimensional sys- 
tems the electrons are localized on a scale £ ~ I. In 
a quasi-one-dimensional geometry however, the localiza- 
tion occurs on a scale £ ~ (kpL±) 2 £ 3> t, where (kpL±) 2 
is the number of transverse channels. Since ( > L we 
expect the electrons to behave diffusivelyO 

In a clean ring we expect the persistent current to be 
on the order of the ballistic current scale, 
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A. The model 

A field-theoretic description of interacting electrons in 
a random potential starts from an expression for the par- 
tition function Z in terms of a functional integral of the 
form 



Z= \ D$,il>]e 



(•5) 



Here ip a (x, r) is a field that depends on the spin la- 
bel a = ±i, position x, and imaginary time r, while 
"0cr(x, r) denotes the conjugate field. The action S[^p,ip] 
includes both the disorder potential and the repulsive 
electron-electron interaction. Since we are dealing with 
a fermionic system, ip and ip are anticommuting (Grass- 
mann) fields; their Fourier transform is defined by 



Vv(x,t) = 




,T) Vv(k, oj n ); (6a) 
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where ui„ = 2tt(ti + |), with n integer, is a fermionic 
Matsubara frequency. The free energy is given by 



= -T(lnZ), 



(7) 



where (• • •) denotes an average over disorder and T is the 
temperature (we use units in which Boltzmann's constant 
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k-Q = 1). To evaluate (lnZ) we use the replica trick, 
which makes use of the identity 
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lnZ= lim — {Z N -1), 



where the partition function 



Z N = 



N 



J{Dir,ip a } exp 
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describes N identical "replicas" of the original system. 
Averaging over a Gaussian disorder distribution intro- 
duces quartic interaction terms between different repli- 
cas; these terms occur in addition to single-replica quar- 
tic terms describing electron-electron interaction. Sub- 
sequently, all these four-fermion terms are decoupled by 
performing Hubbard-Stratonovich transformations. This 
introduces a classical matrix field Q"m with replica labels 
a, (3 and Matsubara frequency labels n, m. The resulting 
action is quadratic in the fcrmionic degrees of freedom tp 
and so that these can be integrated out. The result is 
a classical field theory of the form: 



Z N = J D[Q]e s ^. 



(10) 



The action S[Q] has a continuous symmetry that was 
first recognized for noninteracting electrons by Wegner.Eil 
In the low frequency limit, the Goldstone modes associ- 
ated with this symmetry correspond to particle-hole ex- 
citations ( "diffusons" ) with Q^m corresponding to the 
composite variable - 0"V'm or V'nV'm; an d particle-particle 
excitations ( "cooperons" ) with Q"^ corresponding to 
<<orC^. 

The field Q"^(x) is an element in four-dimensional 
bispinor space, spanned by the charge (particle/hole) and 
spin (up/down) degrees of freedom. A suitable basis in 
this space is formed by the matrices e r <g> Si , where r,i — 
0, . . . , 3. The matrices 
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span the charge space. Note that this definition is dif- 
ferent from the one used in Ref. [2(], since the latter is 
not useful for taking into account the effect of the mag- 
netic flux enclosed by the ring (see Sec. [II C below). This 
modification gives rise to a somewhat different formalism 
at a technical level; we will explicitly indicate differences 
as they occur. For the spin space we use the same repre- 
sentation as in Ref. Eq. It is spanned by 



•so = 
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1 
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S2 = 



1 

-1 



S3 = 







In terms of this basis, Q^m is given by 

3 3 

Qnm = 5Z( er ® Si ) rQnm 

r=0 i=0 



(11) 



The matrix Q can be chosen to be Hermitian: 

Q = Q f . (12) 

In addition, it is required that as an operator in bispinor 
space, Q satisfy the condition 

C T Q T C = Q, (13) 

where C is the charge conjugation operator, given by 

s 2 
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(14) 
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It follows from Eqs. (||) and (|l|) that 
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We expand around the "Fermi-liquid" saddle pointS 

(Qsp)"™ = [(eo + £3) ® s 2 ] S af) S nm sgn(w„). (16) 
An alternative notation for Qsp is 

1 



Qsp = 



-1 



(17) 



where the 2x2 matrix structure refers to the Matsubara 
frequency indices in the following way: 



n, m > 



n < 0, m > 



n > 0, m < 



n, m < 



(18) 



In addition to the slow Goldstone modes there are mas- 
sive modes that can effectively be integrated out by im- 
posing the constraints 



Q 



1: 



tr Q = 0. 



(19) 



These constraints can be eliminated by choosing a suit- 
able parametrization: we will consider fluctuations 



= Q - Qsp, 



of the form 



V 7 ! - OT f - 1 



(20) 



(21) 
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where the unconstrained field q^ n (x) has n > 0, m < 0; 
it can again be expanded in terms of (e r (g> Sj) in exactly 
the same way as Q [see Eq. (|TT|)] . Expanding the action 
S[Q] in terms of the gradient operator and frequency, and 
neglecting terms of order V 4 , uj 2 , and V 2 o-i, we obtain the 
action 



S[Q] = ^ / ( /xlr(V() ( x)r 



2H / dxtr(nQ(x))+5 in t[Q], 



(22) 
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where 



= 5 a p8 nm uj n (e + e 3 ) (g) s - 



The coupling constants G and are defined as 



G 



ttN f D' 



H = 



wNp 



(23) 



(24) 



where Np denotes the saddle point density of states at 
the Fermi energy. Note that GH = l/D. The interaction 
term Si n t[Q] term is defined in Ref. [2(| It is quadratic 
in Q and consists of three parts, with coupling constants 
K s for the singlet channel, K t for the triplet channel, and 
K c for the singlet Cooper channel; the coupling constant 
for the triplet Cooper channel vanishes as a consequence 
of the Fermi exclusion principle. 

Using the parametrization (|2l]) we can write the action 
in terms of an expansion in powers of q as 



S[q] = J2S^[q}. 



(25) 



It should be stressed that although this is an expansion 
in the strength of the disorder, l/kpi, the interaction 
coupling constants K s t c are taken into account to all 
orders for each term in the expansion. 

Higher order gradient-terms can be added to the action 
in a systematic manner .E3 These terms can be shown give 
contributions to the persistent current that are smaller 
than the 0(V 2 ) result by a factor (£/L) 2 , and they can 
consequently be neglected. 

When the behavior of the coupling constants under 
renormalization group (RG) transformations is consid- 
ered, the interaction constant may under certainj-cifc. 
cumstances be strongly renormalized to large values.c3'E£l 
However, since here we are dealing with mesoscopic sys- 
tems, in which the finite system size provides an upper 
bound on the spatial rescaling underlying the RG flow 
equations, these renormalization effects will be unimpor- 
tant for our purposes. 



B. Fluxes and replicas 

The key ingredient in the calculation the nth moment 
of the current distribution, (7™(0)), is the introduction of 
n classes of replicas, each having their own independent 
flux fa. Let the total number of replicas be N — Y^i=i ^i, 
and let the fluxes be assigned as follows: 

a = 1,2, JVi , (iVi + 1), .(JVi+iVfr) , • ■ • , 

(JVi + . . . + N n -i +!),•••, (Ni + . . . + N n ) . (26) 



Y[Z N '(<t>i) = J] [l + ;VilnZ(&) + 0(N?)] , (27) 



we derive that 
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(-T) n nT"q n""o \ JLl jy. 



(28) 



Using the classical path integral representation for (Z N ) 
given in Eq. (|To|) we thus have for the first and second 
moment: 
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(29b) 



C. Fluxes and wave numbers 

Next we turn to the question of which wave num- 
bers are allowed. Since the width of the ring is much 
smaller than its circumference, L± <C L, we neglect exci- 
tations in the transverse direction, i.e. we only consider 
k = (fc,0, 0). At low temperatures in particular, this 
"freezing out" of the transverse degrees of freedom is 
justified since transverse excitations have a much higher 
energy than longitudinal ones. 

We furthermore assume that the vector potential is 
constant over the circumference of the ring. We can then 
'gauge away' the magnetic field, and summarize the effect 
of the magnetic field by means of the twisted boundary 
condition 



tf>{x + L) = e 2 ^l^^){x) 



(30) 



where </>o = 27r/e is the elementary flux quantum, in units 
where % = c = 1. Consequently the allowed wave num- 
bers, with rik an integer, are 
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ip(k) : k = — I n k + -f- 
L \ 0o 
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t/j(k) : k = — [ n k - -f- 
L \ 0o 



(31a) 



(31b) 



Using the identity 



It follows from the definition of e r that for IQnm (p) with 
p = (p, 0, 0) the allowed wave numbers p, with n p an 
integer, are 
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(32) 



Note that for diffusons (r = 0, 3) the wave number p 
depends on (<fi a — <f>@) while in the case of cooperons (r = 
1, 2) we have a dependence on (4> a + <j)p). 



D. The zero-momentum mode 

When a field theory in a finite geometry is consid- 
ered, the zero-momentum mode usually requires special 
treatment .E3 In the case of the nonlinear a model the 
integration over the zero-mode corresponds to an aver- 
age over the group of transformations that leave the ac- 
tion invariant. The moments of the current distribution 
are invariant under this group. However, they are cal- 
culated by means of an expansion around a particular 
saddle point that breaks the symmetry of the group. It 
is for this reason that the zero mode should be omitted 
in all expressions. 

As will become clear later on, the expressions for the 
current moments will in general contain sums over dis- 
crete momenta that are shifted by the flux <j) in the 
case of cooperons, while there is no dependence on <p 
for diffusons. For example, when in Eq. ( |32| ) we choose 
4>a = 4>p — the values of p depend on ip for r = 1, 2 
while they are independent of <fi for r = 0,3. Therefore 
only in the case of diffusons does the zero-mode p = 



occur, and should it be omitted. Later on, in Sec. VD 
we will come back to the question what happens in the 
case of cooperons when the flux <j> is close to a multiple 

of (j) . 



IV. GAUSSIAN THEORY 

A. The Gaussian propagator 

To quadratic order in q the action is given by the Gaus- 
sian part only, 



GLL 2 



r«12(p) j.Mi2,34(p) Uti(P)i 

p r,i 

(33) 



where 1 = {ai,rti}, etc. In the case of diffusons (r = 
0,3), the matrix M is given by 

0,3^^12,34 (p) ^0103 ^a2ft4 ^ril — n2 ,n3 — 714 

X {<Srnn.jP 2 + GH{uJ ni - LU n2 )] 



2-nTGK, 



t}, (34) 



where the interaction term contains the singlet coupling 
constant K s for i — and the triplet coupling constant 
K t for i = 1,2, 3. For cooperons (r = 1, 2) we have 

i,2-Mi2,34(p) = 8 

Ot\ CK3 $0.20-4 +?12 ,713+714 

X {<Wn 3 b 2 + GH(bJ ni - LO n2 )] 



+ S it o S aia2 2ttTGK c }, (35) 



where K c is the cooperon coupling constant. 
The Gaussian propagator is given by 



<r?12(Pl)i^4(P2)> = 



GLL 2 , 



Srs Sij <5 P1 , P2 ^12)34 (Pi)- 

(36) 



Note the slight difference with Ref. |2(j. Another mod- 
ification is that due to the symmetry property ( |l5|) we 
have 

(r9l2(Pl)3- S 934(P2)) = (r<7l2 (Pi )i<?34 (P2)) 

= (I-A2 (Pi) J s <&(P2)> (37) 

The inverse M" 1 of the matrix M has been calculated 
in Ref. We cite the result, in terms of the diffusive 
propagators 

£>»(p)= [p 2 + GHn n ]~\ (38) 
V s n t (p)=[p 2 + G(H + K s ^n n ]-\ 
&Vf{p)=Vf{p)-V n {p). (39) 



For diffusons we have 

0,3-^12,34 (p) = ^a 1 a 3 ^a 2 a 4 ^n 1 ~n 2 ,n 3 —n 4 



x S S nirl3 2? ni _„ 2 (p) 



AVI 



(P) , (40) 



ni — n% 
while for cooperons 

1,2-^12,34 (P) = ^aia 3 <5a 2 c<4 ^111+112,113+114 

X — 712 (p) L) ni — n 2 (p) ^n 3 — 714 

Si,0 S aia2 2-kTGKc I 



where 



1 + 2nTGK c f ni+n2 (p) 



/m+n 2 (p) — S n ^ ni+n2 V ni+n2 (p). 



(P) 

(41) 
(42) 



Note that contributions related to electron-electron in- 
teractions (i.e., proportional to if s ,t,c) require that all 
replicas be the same: a\ = 012 = &3 — 
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B. Logarithmic suppression of K c 



where the allowed wave numbers are 



The interaction part of the cooperon propagator (i.e., 
the part proportional to K c ) is suppressed by a factor [1 + 
2irTGK c fn 1+n2 (p)]~ 1 . It is important to estimate the 
quantitative effect of this term. The term f n (jp) diverges 
logarithmically, as 



2nTGK c f n (p) = 8ttTN f V 



f 



m=0 
,1/t 

2N F V / du 
Jt 



Dp 2 + 2-nT{n + 2m) 



1 



' Dp 2 + Lu n + 2uj t< 



2N F Vhi 



Dp 2 + u> r , 



Dp 2 +u) n + T 



(43) 



The upper and lower bounds on the frequency integral 
are provided by the average collision frequency and the 
temperature, respectively. 

Since the electron-electron interactions are screened 
and therefore short ranged, they can be approximated 
by a point interaction with strength V . From the def- 
inition of K c in Ref. || it follows that K c = ttN^V/2, 
or GK C = 2N F V/D. A typical valueO for the electron- 
electron interaction is N F V ~ 0.3. 

The dominant contribution to the persistent current 
comes from wave numbers and frequencies on diffusive 
scales, p ~ 2-k/L and w„ ~ D(2tt/L) 2 , and T < T 
with T = D/L 2 , so that the argument of the loga- 
rithm has magnitude (L/2tt) 2 /Dt = 3(L / £) 2 / (2n) 2 > 1. 
Using N F V ~ 0.3 and Ljt ~ 140 we estimate-, that 
{l + 2irTGK c f n (p)] ~ 5, in agreement with Eckern,E3 who 
obtained the same result by summing an infinite series of 
Feynman diagrams. 

In the remainder of the paper we will for each occur- 
rence of the cooperon propagator replace N F V by its ef- 
fective value N F Vcs ^ 0.2N F V ~ 0.06, and drop the 
denominator [1 + 2nTGK c f(p)\. 

It is worth noting the close connection between the log- 
arithmic suppression described here and the logarithm 
appearing in BCS theory. Both result from the summa- 
tion of an infinite series of diagrams involving multiple 
interactions. In the case of a phonon-mediated attractive 
electron-electron interaction, there is an enhancement 
rather than a suppression of K c , which for low enough 
temperatures leads to destabilization of the Fermi liquid 
ground state and the emergence of superconductivity. 



C. Calculation of </) (2) 

To calculate (I)® in Gaussian approximation we note 
that 



dS® [q] -2 



d(j) GLL\ 



EEE( 2p ) Jr UMp) Wi2(p), 

p 12 ri ™ 

(44) 
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(45) 



All N = Ni replicas have flux <j). Due to presence of the 
factor (dp/d(j>) only cooperons (r = 1,2) contribute to 
the average current. The result is 



T 
2 



p 12 r=l,2 



dp 



(46) 



In the replica limit (N — > 0) only terms proportional to N 
survive. It follows that only the interacting part of M _1 
contributes, since it is diagonal in the replica labels; the 
nonintcracting part of M _1 yields a term proportional to 
TV 2 and consequently vanishes in the replica limit. The 
previous expression for the average current can be written 
as 



with an average free energy given by 



<n(</>))( 2 ) = 2n f vtY,J2 



p uj>Q 



Dp 2 +u 



(47) 



(48) 



where u> = 2nTn, with integer n > 1, and p = (2tt/L) 
x(n p — 2<j)/(f>o). This is exactly the Hartree-Fock result 
obtained in Ref. [ll], and the corresponding average cur- 
rent at T = is given by 

m) * * i D ( y ^(w/0o) (49) 

V 37r / ^-^ m 6 

where Id is the diffusive current scale defined in Eq. (^) 
m Sec. 0. This result was obtained by exploiting the 
fact that (I((f>)) is a periodic function of <f> with period 
4>a /2. In the calculation of the Fourier components the 
sum over wave numbers gets replaced by an integral (see 
Ref. [ll]). An alternative calculation that does not rely on 
Fourier transformation of the <p dependence shows that 



E 

m— 1 



irr 



dss 



(sinhs)(sin 4-7T0) 
(coshs — cos47T0) 2 

(50) 



Note that the integral exists for all cf>. 

When the "bare" interaction V rather than the ef- 
fective interaction V c h is used, the experimental and 
theoreticalliil values for the average current (I) are in rea- 
sonable agreement. The loga rithm ic suppression of the 
interaction discussed in Sec. IV B severely worsens the 
agreement between theory and experiment, leading to a 
discrepancy by a factor of 5. 
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It is however perfectly consistent to keep the suppres- 
sion term. As was already noted earlier, the nonlinear 
a model provides an expansion in terms of the disorder 
strength, 1/kpl. The inversion of the asymmetric matrix 
M performed in the calculation of the Gaussian propa- 
gator, however, amounts to an infinite resummation in 
terms of the interaction constant K c . Therefore in a cer- 
tain sense the electron-electron interaction is taken into 
account exactly. 



D. Calculation of (J 2 ) (2) 

The second moment of the current distribution can be 
written as a sum of a disconnected and a connected part, 



(I 2 ) = (I) 2 + (I 2 ) C 



(51) 



According to Eq. ( 29b ) there are two contributions to the 
second moment (7 2 )^ . We first consider the term 



d 2 S®[q] 



-4 
GLL\ 



Z^Z^Z^ /m, f)A„ r 



p 12 r,i 



I J 12 (p)* 

(52) 



Here the wave numbers are given by 



2tt 

P = T 



V 



</>0 



Va-2 

4>o 



(53) 



so that (dp/d4>i)(dp/d(f>2) can only be nonvanishing if 
either (i) <fi ai — <pi and <p a2 = 4>2, or (ii) 4> a2 = 4>\ and 
4>ai = 4>2- In both cases we necessarily have a\ ^ a.2, 
so that only the noninteracting part of M _1 can con- 
tribute. Furthermore it can easily be verified that in 
both cases the summation over replica indices yields pre- 
cisely the factor N1N2 that is needed for the term to sur- 
vive in the replica limit. Note that (dp/d(f>i)(dp/d(j>2) = 
±e 2 (2ir/L) 2 , with a plus sign for diffusons (r = 0, 3) and 
a minus sign for cooperons (r = 1,2). 

A second contribution to (J 2 )® is of the type 



\GLL\ 



E v \ - dpi dp! 
z^ z^ f)^, pi^ 



/ dS®[q}dS®{g} \ (2) 

\ 9(1)1 d ^ 2 I x/ pl , P2 i 2 34r s ,y 

X (2 Pl )(2p 2 ) (Ul2(Pl)U*12(Pl)MP2)lqUP2)f ■ 

(54) 

When applying Wick's theorem to the four-point corre- 
lation function (qq*qq*)® there are two ways of pairing 
q with q*: 



{ l r qi2(pi)l.q*12(Pl)l<l3i{P2)iq34(P2)) 

(pi) i-M^M 

2 



(2) 



12,12 U'l 

-1 



M- 



12,:S4 



{pi)] 



(55) 



Comparison with Eq. ( [iq) shows that the first term on 
the right hand side corresponds to the disconnected part 
((/) (2) ) 2 . The second ter m on the right hand side how- 
ever yields a contribution to the connected part (J 2 )? . 
As M^2 34(p) ~ 5 aia3 5 a2ai it is again only the nonin- 
teracting part of 34(2^) that can give rise to a term 
proportional to NiN 2 . 

Again we can identify our result with a_ result ob- 
tained using standard perturbation theory.E3 Diffusons 
and cooperons can be shown to give identical contribu- 
tions. It can be verified that we can write 



d 2 M</W)> 



with 



WW)) = 



AT 



EE ln 

p u>Q 



Dp 2 



(56) 



(57) 



where to — 2nTn, with integer n > 1, and p = (2n/L) 
x (n p + (<f> + <j)')/(j) ). This is exactly the result obtained 
in Ref. [LOL which at T = is given by 



W) 2 



D 



'2\/2 N 
~3V 





'2vr0" 


sin 









(58) 



We stress again that this result is independent of the 
electron-electron interaction. Numerically, it falls short 
of explaining the experiment of Ref. |i| by one or two 
orders of magnitude. 



E. Persistent current distribution is Gaussian 

In this section we will show that the higher moments of 
the current distribution, (J™)® with n > 3, are exactly 
those of a Gaussian distribution, provided that we use 
only the Gaussian part, of the action when evalu- 

ating them, i.e. we approximate (• • •) by (• • -)^ 2 \ In the 
remainder of this section we will suppress the Gaussian 
superscript for notational convenience. 

Let us start by observing that 



S[q]_ 



d" 



oS[q] 



Til .712 >0 



In this expression the fluxes <fi* denote a permutation of 
the fluxes {</>i}™ =1 - The right hand side leads to corre- 
lation functions of the type ((qq* ) ni (qq*) n2 ), which are 
factorized using Wick's theorem. 
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To facilitate the combinatorics of the problem it is 
helpful to introduce a symbolic notation. For the average 
current we have 



with 



(-0 



{qq 



(qq*). 



(60) 



The symbol corresponds to a differentiation with 

respect to the flux, dS/d<f>, which produces a term qq*; 
furthermore this symbol implies a term (2p)(dp/d(j>) (see 
Eq. ([46])). The line under qq* denotes a pairing of q and 
q* that produces a propagator M~ l . 

For the second moment we have in symbolic notation 



(I 2 ) = (qq*qq*) + (qq*) 

= (qq*) (qq*) + (qq*qq*) + (qq* 



(61) 



On the last line of this equation, the first term corre- 
sponds to the disconnected part of (I 2 ) which is just the 
square of the average current; the remaining two terms 
constitute the connected part (I 2 ) c , where the second 
term corresponds to (dS / d<pi)(dS / dfo) and the third 
term corresponds to (d 2 S/d(j)id(j)2)- Both contributions 
to (I 2 ) c contain a factor (dpi / d4>i)(dp2 / dfa) ■ The first 
two terms arise from the two different ways in which q 
can be paired wit h q* i n a product of the type qq*qq*, as 



discussed in Sec. IV D 



To understand what happens for the moments n > 3 it 
is important to recall that any wave number p can depend 
on only two different fluxes (see Sec. Ill C ). When two q's 
are paired, their replica indices must be the same since 
M q 1 1 q 2 q 3Q4 ~ S aia2 5 a3a4 . Consequently, when terms 
containing three or more fluxes are paired into a con- 
nected cluster, the corresponding expression will contain 
a term, 



dp dp dp 



= 0, 



(62) 



which vanishes by definition, since p cannot depend on 
three fluxes. In terms of the symbolic notation used here 
we have, e.g., 



(qq*qq*) = (qq*qq*qq*) = 0. 



(63) 



It follows that when Wick's theorem is applied to a gen- 
eral term of the type 



{qq* ■ ■ ■ qq* qq*- ■ ■ qq* 



(64) 



it factorizes into a product of connected terms that de- 
pend on either one or two fluxes. There are three such 
terms, and all of them occur in the expressions for (I) 
and (I 2 ) c - Working out the combinatorics of the possible 
ways of pairing we obtain 



(I' n )= E C km ((£*)) 

k,m>0 
2k+2m<n 

(( qq*qq* )) ((qq*)) 



n— 2k— 2m 



(65) 



■'km 



k\ m\ (n — 2k — 2m) 



(5 



1 \ k-\-m 



(66) 



The characteristic function p(z) is defined in terms of the 
persistent current distribution function p(I) as 



p(z) 



Substituting Eqs. (|6J]) and 
sion we obtain 



e ^) = E^ 

1 n=0 



p(z) 



exp 
exp 



z{qq 
z(I) 



(67) 

for (I™) into this expres- 

!) + h z2 (( qq*qq* ) + (qq3 . 



W(I 2 



(68) 



From this expression we read off that, as expected, the 
first two cumulants of the distribution are (7) and (I 2 ) c , 
respectively, while all higher cumulants vanish. This is 
the signature of the Gaussian distribution. 



V. BEYOND GAUSSIAN APPROXIMATION 

As we saw in the previous section, the only part of the 
action that depends on the flux </> is the gradient term 



dx tr(VQ(x)) 2 



(69) 



In evaluating averages of the type (qq* ■ ■ ■ qq*) we used 
the full Gaussian part S® [Q] of the action, i.e., all terms 
of order 0(q 2 ) including the interaction terms. This is 
the reason why the Gaussian result for (I) can depend on 
the interaction coupling constant K, even if Sy [Q] does 
not. 

In what follows we will take into account terms of 
higher order in the expansion in powers of q of SV[Q]- 
Averages (qq* ■ ■ ■ qq*) however will still be evaluated in 
Gaussian approximation with a weig ht factor exp(S<%]), 
so that Wick's theorem can be used to express the aver- 
age in terms of products of Gaussian propagators, as was 
done in the previous section in the case of (I 2 ). Since 
the disorder is relatively weak, it is expected that using 
the full average (• • •) instead of the Gaussian approxima- 
tion (■ • -)( 2 ) would simply lead to small corrections to the 
terms computed here using the Gaussian average. 



A. Calculation of (J) (4) 

The first term in the expansion beyond Gaussian order 
in the expansion of SV [Q] in powers of q is given by 



S v [q] - 2G(LL'i) 



Wee e e Ep 2 



(70) 



p fci,k 2 1234 {r a } {i a } 
X l\qi 2 (p + ki) \\ q* 2 (hi ) % (734 (-P + k 2 ) l r \qU(k2) 

x tr(e n 4 2 e r 3 4 4 ) tr ( s n 4 2 S i3 s jj- ( 71 ) 
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It should be noted that there is a very similar term 
in which the role of q and is interchanged. It can 
be verified that this term yields an identical contribu- 
tion, so that it is sufficient to multiply the final result 
by a factor 2. The quadruplets {rir 2 r^r^) for which 



tr(e ri e 



4,) = i 



arc 



r 2 Cr 3 c r 4 1 

(0000) (1111) (2222) (3333) 

(0220) (1331) (3113) (2002) 

(0011) (2233) (3322) (1100) 

(0231) (1320) (2013) (3102) 

while the trace vanishes for all other quadruplets. Note 
that the allowed quadruplets all have the following prop- 
erty: 



and 



















+ 




+ 


) 


+ 




+ 


> (72) 








r 4 




T 2 


. + _ 




+ 




+ 


? 


+ 




+ 


■ (73) 



















The contribution of [q] to the average current (J) 
reads 



T 



/ dS®[q] 
\ 90 



r EEEEE( 2 ?) 

p fei.fca 1234 {r a } {i a } 



dp 



x (Ji?ia(p + *i)^952(*i)^te4(-P + Aa)ji^4(*a)> > 

(74) 

with the allowed values for p given in Eq. ( |45| ) . As in the 
preceding section, there are two ways of pairing q with 
q* . Pairing q\2 (p + fci) with q^ 2 {ki) yields a Gaussian 
propagator proportional to S Pt o, which together with the 
factor 2p leads to a vanishing contribution. The alterna- 
tive pairing of qi2 (p + fci) and q$ 2 (k 2 ) gives rise to the 
following expression: 



(I) 



(4) 



2T EEEE 

fci ,k 2 1234 n ,r 2 ii ,i 2 



(k 2 - fci) 



dp 



(75) 



Note that r± — r 4 , r 2 = r 3 , a 2 = 0.4,, and n 2 = TI4. The 
only contribution that survives in the replica limit is the 
one where all four replica labels are the same, which is 
only possible if for both propagators M^ 1 we take into 
account the interacting part only, as this part is com- 
pletely diagonal in the replica labels. We then have 



w (4) - E 



(76) 



To proceed it is necessary to analyze the dependence 
of fci, k 2 , and p on the flux <fi and the indices r\ and 
r 2 ■ In the Gaussian case we saw that only cooperons can 
contribute; here however the situation is slightly more 
complicated. We have 



2tt 

T 



fci 





( 












n p + 




+ 




+ 








\ 




ri 




J 





(77a) 



2tt 

T 



( 


■ - 


\ 






2(j) 










+ 




V 


. . 


r 2 I 



(77b) 





f 


- " 


\ 


2tt 








T 




+ 








. . 


ri J 



(77c) 



From these relations it can be seen that the allowed 
quadruplets of the form {r\r 2 r 2 r\) fall into three classes: 
(i) dp/d(t) = -2 for (nr 2 ) = (02) or (13); (ii) dp/dcj) = +2 
for (nrj) = (20) or (31); (iii) dp/d<j> = for r x = r 2 . It 
can be seen that the pairs (02), (20), (13), and (31) all 
give identical contributions to (/)". For example, the 
contributions of (02) and (13) are equal since the ex- 
pression for (J)® is invariant under the transformation: 
fci — > -k 2 , k 2 — > —fci. 



We are now ready to evaluate Eq. ( 76 ) . Focusing on the 
case (rir 2 ) = (02) and introducing a rescaled frequency 
y = (u>/D)(L/2tt) 2 we obtain 



(78) 



We have written the prefactor of -F(^) in a way that 
allows for an easy comparison with the Gaussian order 
result, but it should be noted that this expression does 
not depend on the strength of the disorder, as can be seen 
by using D — vp£/3. The flux dependence is contained 
in the dimcnsionless expression, 



F{4>) = dy 



n 2 + 20/ <f> 



<>2 

2tt / ds 
Jo 



(™ 2 + W/M 2 + V ^ n\ + y 

sin47T(/) / 1\ 

— v - , coth s - - , (79) 

cosh 2s — cos47T0 V s J 



where we have used a change of variable to s = TT^/y. 
The zero-momentum term rt\ = h as bee n excluded, the 
rationale for which was given in Sec. Ill D . The integral is 
convergent for all values of 4>; this would not be the case 
if the zero-mode contribution had not been subtracted, 
since then there would be a small-frequency divergence. 
We will now argue that the term (I)^ and all other 
higher-order terms, are small compared to the Gaussian 
contribution (1)^ . 
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B. Magnitude of general higher-order terms for (I) 

In the preceding section we have performed the cal- 
culation of the first contribution to the average current 
beyond the Gaussian order. We will now generalize this 
calculation to estimate the magnitude of contributions of 
arbitrary order to (I) and (I 2 )- 

Our aim is to write all expressions as a product of a 
dimensionful prefactor, i.e., a combination of the param- 
eters D, L, Lj_, Np, and V , and a dimensionless sum over 
rescaled wave numbers and frequencies. Consider there- 
fore the following contribution to the average current: 



(2r, 



= -T 



\ d<t> 



(80) 



First, there is a factor l/[G(LL 2 L ) 2m 1 ], associated with 
S&^lq], Next, the terms (2p) and (dp/ dip) each yield 
a factor (2tt/L). The differentiation of S^[q] yields 
a (2m)-point correlation function in terms of q; apply- 
ing Wick's theorem we obtain a product of m Gaussian 
propagators M _1 . Each propagator has a factor GLL± 
associated with it. 

The propagators can either contribute their noninter- 
acting or their interacting part; the latter is diagonal in 
the replica labels. The number of interacting contribu- 
tions is determined by the requirement that due to the 
replica limit there should only be a single independent 
replica index. This can only be achieved by choosing the 
interacting part of all propagators, which is proportional 
to [D(k,oj)] 2 . There is a factor (L/27r) 4 coming from a 
rescaling of [D(k, uj)] 2 , and a factor GK = N F V/D repre- 
senting the interaction. For simplicity, wc have neglected 
the logar ithmic correction term that was discussed in 
Sec. [V C . Collecting all factors we find that 



(I) 



(2r, 



e(GLLl) m (T )P 
G(LLi) 2 ™- 1 



(81) 



The factor (Tq) p , with To = D/L 2 , comes from a rescal- 
ing of the frequency variables. Rather than by working 
out the structure of the frequency summations in detail, 
we can determine p by dimensional analysis. From the 
requirement that (J) be proportional to a frequency, it 
follows that p = m + 1, so that 



L 

1~ 



1 



NpDL 



(82) 



Note that m = 1 corresponds to the Gaussian result of 



Sec. [VC 



To determine the value of the dimensionless combi- 
nation (L I L j_) 2 / (N F D L) we note that in three dimen- 
sions the density of states at the Fermi energy is given 
by N F = mk F /-K 2 , while D = vpl/Z, so that N F DL = 



(k F £)(k F L)/3ir 2 ~ 10 6 . This is so large that in spite of 
the large factor (L/L±) 2 ~ 10 4 we have 



i 



N F DL 



10" 



(83) 



Thus we conclude that the higher-order terms are small 
compared to the Gaussian result. 

Although contributions thus seem to be smaller by a 
factor of about 10~ 2 for each higher order, a simulta- 
neous rapid growth occurs of the combinatorial factors 
arising, e.g., from the different ways of pairing q's when 
applying Wick's theorem. The perturbation expansion in 
terms of the disorder therefore at best provides us with 
an asymptotic series, a situation not unfamiliar in other 
contexts. 



C. Magnitude of general higher-order terms for (I 2 ) 

The calculation of the magnitude of higher-order con- 
tributions to the second moment (I 2 ) proceeds in a simi- 
lar fashion. As we have seen above in Sec. [VD there are 
two types of contributions. First we have 



(T 2U2 m)=T 2 j d 2 S^[q] 

K 0000' 



(84) 



A difference with the calculation of (J) 2ra is that the 
replica limit now requires the existence of two indepen- 
dent replica indices. This can be achieved by choosing 
the noninteracting part of one of the propagators, while 
using the interacting part of all (m — 1) others. The 
noninteracting propagator contributes a factor (L/2w) 2 
instead of (GK )(L/2tt) 4 , so that 



(I 



2\ (2n 



e 2 (GLL 2 ± ) m (T ) p 



G(LL 



2 \2m-l 



(85) 



In this case the requirement that (I 2 ) be proportional to 
a frequency squared, again yields p = m + 1, so that 



(I 



2\ (2m) 



e^D 2 



(N F V) ( — 



1 



-i m—1 



N F DL 



(86) 



The second type of contribution to (I 2 ) is 

•dS^[q] 0S (2m2 >[q] 



00 



e 2 D 2 
L A 



(N F V) mi 



L 



1 



L± J NpDL 



m 1 +m,2 — 2 



(87) 



In this case the Gaussian result corresponds to m — 1 and 
mi = mi = 1, respectively. Again we conclude that the 
corrections for (J 2 ) are small compared to the Gaussian 
result. 
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D. Convergence of the dimensionless sums 



VI. DISCUSSION 




In the preceding subsections we argued that all terms 
beyond Gaussian order are small compared to the Gaus- 
sian result. It is essential to verify that the corresponding 
dimensionless sums over wave numbers and frequencies 
are convergent. Based on the structure of the calculation 
of (I)' 4 ' we expect that any higher-order term (/™)( 2m ) 
will contain a dimensionless, flux-dependent part given 
by an expression of the type 



Here x denotes some linear combination of dimensionless 
fluxes <fii/4>o, an d the integers a and b, satisfying a + 
b = to, are determined by the detailed structure of the 
allowed sets {r^} of particle-hole indices. In the case of 
(iy 4 ) we have x = 2(j>/^>o and a = b = 1. The integral 
over y is always convergent for y — > oo; for generic values 
of x there is also convergence for y — > 0. Only when x is 
close to an integer a divergence may occur. Let us restrict 
ourselves without loss of generality to values \x\ <C 1. It 
is easy to show by a rescaling of y that F%(x) <~ x 2 ~ a as 
x — > 0. In other words, only when a = 1 does the moment 
of the persistent current distribution under consideration 
go to zero in the limit of vanishing enclosed flux. On the 
other hand, it is clear from symmetry considerations that 
F%(x) vanishes whenever a; is a multiple of h. 

In principle it is possible that the case a < 2 never 
occurs, even at high orders in perturbation theory, but 
there is another reason why the apparent discontinuity 
at x = should not be too alarming even if it would 
occur. To see this, one should realize that the electrons 
are localized on a spatial scale £ = £(kFL±) 2 ^> L. Val- 
ues |a;| <C 1 correspond to length scales L/x L that 
will be on the order of £ if x ~ (L/£)(kpL±)~ 2 . On 
this scale there no longer exists diffusive behavior. In 
the context of the nonlinear a model diffusive behav- 
ior corresponds to an ordered state breaking the sym- 
metry of the model. Goldstone modes and their associ- 
ated diffusive propagators only occur in the ordered state. 
This means that for exactly those values of x for which 
Fq(x) diverges, the nonlinear a model description of the 
quantum-mechanical problem breaks down due to local- 
ization effects; the diffusive poles that are the source of 
the divergence of F\[x) are replaced by correlation func- 
tions that are exponentially damped in space, and the 
divergence disappears. It should be noted that the above 
conceptual considerations are somewhat academic, in the 
sense that in the actual experimentsB □ the temperature 
provides a natural cutoff for all possible divergences. 



In this paper we have shown how a nonlinear a model 
approach to the persistent current problem provides us 
with an efficient and unambiguous method for calculating 
moments of the persistent current distribution for a dis- 
ordered ring. We thus have developed a valuable new tool 
for studying the combined effect of quenched impurities 
and electron-electron interaction on the magnitude of the 
persistent current at high orders in perturbation theory. 
More traditional diagrammatic approaches so far have 
failed to provide a reliable method to analyze the conse- 
quences of high order interactions between slow modes. 

Our main result is that higher-order terms order are 
small compared to the lowest order Gaussian result (see 
Sec. |y|). Furthermore, we showed that the Gaussian ap- 
proximation within the nonlinear a model reproduces the 
diagrammatic results of Refs. [Tc|-|l2]. It follows that the 
experiments of Refs. ||J^ cannot be explained purely in 
terms of interacting Goldstone modes, i.e., diffusons and 
cooperons, but that an additional theoretical ingredient 
is needed. 

Using only the dominant Gaussian part of the action 
to calculate the higher moments of the persistent current 
distribution, we found that this distribution is Gaussian 



(see Sec. IV E). This feature is not expected to survive 
in a more complete theory for the tail of the distribution. 

The nonlinear a model is essentially a small wave- 
number and small frequency expansion. However, it is 
known that "nonhydrodynamic" contributions may be 
important, an example being the density [Expansion of 
transport coefficients for quantum systems.E2l It is there- 
fore desirable to investigate the validity of the "hydrody- 
namic" approximation, by rederiving some of the known 
results without performing a small wave number expan- 
sion. The relative importance of contributions coming 
from wave numbers in the diffusive (k ~ 1/L) and ballis- 
tic (k ~ l/l) range can then be determined. An anal- 
ysis along these lines will be presented in a separate 
publication.!! 1 ! 

We believe that the method presented here can be 
used with some confidence, as in the simplest possi- 
ble (Gaussian) approximation we exactly reproduce well- 
known and widely accepted results obtained using stan- 
dard diagrammatic perturbation theory; this was shown 
in Sec. IV. To illustrate the efficiency of our approach we 
note that the derivation of the equivalent of our Gaus- 
sian results using diagrammatic perturbation theory re- 
quires three separate infinite resummations: (i) in stan- 
dard Born approximation the dressed Green's function 
contains an infinite number of impurity lines; (ii) the dif- 
fusons and cooperons used as building blocks in the dia- 
grammatic approach involve infinite ladder summations; 
(ii) the logarithmic suppression of the electron-electron 
interaction in standard perturbation theory is obtained 
by means of an infinite ladder summation of cooperons 
alternated with e-e interaction vertices.EHl 
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An interesting observation is that, as opposed to the 
higher-order contributions, the Gaussian results are in- 
dependent of the width L± of the ring, or alternatively 
the number of channels (kpLx) 2 . Although detailed cal- 
culations rapidly become more complicated beyond the 
Gaussian approximation, it is relatively easy to give a 
naive estimate of the magnitude of higher-order terms 
(see Sec. 0). It was tentatively concluded in Sec. [v] that 
the higher-order terms are small compared to the Gaus- 
sian results. 

All calculations in this paper were performed using the 
grand canonical ensemble, in which the particle number 
fluctuates around an average value determined by the 
chemical potential. Using the method of Altshuler et 
al.fA it would be in principle be possible to calculate 
corrections to the grand canonical results due to a fixed 
number of electrons in the ring. In particular at higher 
order in perturbation theory, it would be interesting to 
assess whether these corrections are negligible or not. 

Our results indicate that although the Gaussian result 
is independent of the thickness L± of the ring, the higher- 
order terms become more important as L± is decreased, 
if all other parameters are kept fixed. At the same time 
however, decreasing L± brings the system closer to a one- 
dimensional system, in which electrons are localized and 
the description in terms of diffusive modes breaks down. 
Thus it is hard to estimate the dependence of the persis- 
tent current on L±, and it would be desirable to experi- 
mentally investigate it. 



ACKNOWLEDGMENTS 

We would like to thank Dietrich Belitz, Bob Dorf- 
man, Manher Jariwala, Arnulf Latz, Raj Mohanty, and 
Masao Yoshimura for stimulating discussions and com- 
ments. This research was performed under NSF Grant 
No. DMR-96-32978. 



1 F. Hund, Ann. Phys. (Leipzig) 32, 102 (1938); 5, 1 (1996). 

2 D. Mailly, C. Chapelier, and A. Benoit, Phys. Rev. Lett. 
70, 2020 (1993). 

3 L. P. Levy, G. Dolan, J. Dunsmuir, and H. Bouchiat, Phys. 
Rev. Lett. 64, 2074 (1990). 

4 V. Chandrasekhar, R. A. Webb, M. J. Brady, M. B. 

Ketchen, W. J. Gallagher, and A. Kleinasser, Phys. Rev. 

Lett. 67, 2578 (1991). 
5 P. Mohanty, E. M. Q. Jariwala, M. B. Ketchen, and R. A. 

Webb, in Quantum Coherence and Decoherence, edited by 

K. Fujikawa and Y. A. Ono (Elsevier, New York, 1996). 
6 N. Byers and C. N. Yang, Phys. Rev. Lett. 7, 46 (1961); 

M. Biittiker, Y. Imry, and R. Landauer, Phys. Lett. 96A, 



365 (1983); R. Landauer and M. Biittiker, Phys. Rev. Lett. 
54, 2049 (1985). 

7 R. Berkovits and Y. Avishai, Europhys. Lett. 29, 475 
(1995); G. Bouzerar and D. Poilblanc, Phys. Rev. B 52, 
10772 (1995); G. Chiappe, J. A. Verges, and E. Louis, Solid 
State Commun. 99, 717 (1996). 

8 H. Kato and Y. Yoshioka, Physica B 212, 251 (1995). 

9 H. Bouchiat and G. Montambaux, J. Phys. (France) 50, 
2695 (1989). 

10 H. F. Cheung, E. K. Riedel, and Y. Gefen, Phys. Rev. Lett. 
62, 587 (1989). 

11 V. Ambegaokar and U. Eckern, Phys. Rev. Lett. 65, 381 

(1991) . 

12 U. Eckern, Z. Phys. B 82, 393 (1991). 

13 U. Eckern and A. Schmid, Europhys. Lett. 18, 457 (1992). 

14 R. A. Smith and V. Ambegaokar, Europhys. Lett. 20, 161 

(1992) . 

15 M. T. Beal-Monod and G. Montambaux, Phys. Rev. B 46, 
7182 (1992). 

16 P. Kopietz, Phys. Rev. Lett. 70, 3123 (1993). 

17 G. Vignale, Phys. Rev. B 50, 7668 (1994). 

18 A. Miiller-Groeling and H. Weidenmiiller Phys. Rev. B 49, 
4752 (1994). 

19 A. Schmid, Phys. Rev. Lett. 66, 80 (1991). 

20 F. von Oppen and E. Riedel, Phys. Rev. Lett. 66, 84 
(1991). 

21 B. L. Altshuler, Y. Gefen, and Y. Imry, Phys. Rev. Lett. 
66, 88 (1991). 

22 A. A. Abrikosov, L. P. Gorkov, and I. E. Dzyaloshinsky, 
Methods of Quantum Field Theory in Statistical Physics 
(Dover, New York, 1975). 



_•:{ 



R. A. Smith, Ph. D. thesis, Cornell University, 1993. 



24 F. Wegner, Z. Phys. 35, 207 (1979). 

25 A. M. Finkel'stein, Zh. Eksp. Teor. Fiz. 84, 168 (1983) 
[Sov. Phys. JETP 57, 97 (1983)]. 

2e D. Belitz and T. R. Kirkpatrick, Rev. Mod. Phys. 66, 261 
(1994). 

27 D. J. Thouless, Phys. Rev. Lett. 39, 1167 (1977). 

28 B. L. Altshuler, V. E. Kravtsov, and I. V. Lerner, in Meso- 
copic Phenomena in Solids, edited by B. L. Altshuler, P. 
A. Lee, and R. A. Webb (Elsevier, Amsterdam, 1991). 

29 J. Zinn-Justin, Quantum Field Theory and Critical Phe- 
nomena (Clarendon, Oxford, 1989). 

30 T. R. Kirkpatrick and J. R. Dorfman, Phys. Rev. A 22, 
1022 (1983); T. R. Kirkpatrick and D. Belitz, Phys. Rev. 
B 34, 2168 (1986). 

31 H. J. Bussemaker and T. R. Kirkpatrick (unpublished). 



12 



